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I derive the space-time distribution of energy and momentum deposited by a fast parton traversing 
a perturbative quark-gluon plasma by considering the fast parton as the source of an external color 
field interacting with the medium. I include the medium's response to screen the fast parton's color 
field by incorporating dielectric functions and compare to the unscreened result. 
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I. INTRODUCTION 

Quantum chromodynamics (QCD) predicts a transi- 
tion in sufRciently hot and/or dense nuclear matter from 
colorless hadrons to a state of deconfined quarks and glu- 
ons known as the quark-gluon plasma (QGP) [H- It is 
believed that this transition has been observed in ex- 
periments performed at the Relativistic Heavy-Ion Col- 
lider (RHIC) (see, for example, [1]). Fast partons (par- 
tons with velocity approaching the speed of light) created 
by hard transverse scattering in the early moments of a 
heavy-ion collision are a useful probe in understanding 
QGP physics. Of particular importance is the process 
of jet quenching in which fast partons lose energy and 
momentum by interacting with the surrounding medium 
(see, for instance, I, i, i, S 0, S i, ) • 

While much effort has been devoted to understanding 
the dynamics of jet quenching, the question of how the 
energy and momentum deposited by the fast parton af- 
fects the bulk behavior of an evolvin g Q GP is relatively 
new (see, for example, [11 [11 [11, QIll [11, [13, EH)- 
One promising approach is to treat the fast parton as the 
source of energy and momentum coupled to the hydrody- 
namic equations of the medium. There is strong evidence 
[20I . [2l| that the matter created in heavy-ion collisions at 
RHIC can be described by nearly ideal hydrodynamics 
and it was argued in [Tsj that if the energy and momen- 
tum deposited by a fast parton thermalize quickly the 
resulting disturbance should evolve hydrodynamically as 
well. 

Casalderrey-Solana et al. [ij] suggested that since the 
fast parton travels through the medium at supersonic 
speed, one might expect a Mach cone to develop in the 
surrounding hydrodynamic medium. They proceeded to 
show that if one uses an appropriately chosen supersonic 
sound source in a linearized hydrodynamic evolution one 
indeed observes a propagating Mach cone. In their work, 
the explicit calculation of the source term generated by 
a supersonic parton was not addressed. In a later work 
I13I Chaudhuri and Heinz examined a two-dimensional 
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hydrodynamics simulation of an expanding QGP with a 
phenomenologically motivated ansatz as their fast parton 
source term. Here the authors failed to find the formation 
of a Mach cone except for seemingly unphysical values of 
energy loss. Much recent work [ll, [ll| has been focused 
on examining a supersonic color charge passing through 
an A/" = 4 supersymmetric Yang-Mills plasma (in both 
of the above citations the authors find a Mach cone in 
the resulting dynamics). Such efforts are motivated by 
the idea that the system is a useful model of a real QCD 
plasma. 

As the above examples illustrate, there have been sev- 
eral insightful attempts at understanding the influence of 
fast partons on the bulk evolution of an expanding QGP. 
However, the explicit evaluation of the distribution of 
energy and momentum deposited by a supersonic parton 
traversing a QCD plasma has remained an open ques- 
tion. It has been pointed out [ll|, [l^l that the medium's 
response to a fast parton is sensitive to the specific form 
of energy and momentum deposition. For example, the 
hydrodynamic source vector of a parton moving at con- 
stant velocity, u, can be written in Fourier space as (Toj 

J(cj,q)= [vL(j)u{uj, q^) +iq(t)q{uj,q^)\2TT6{uj -u-ci). (1) 

The medium's response then becomes sensitive to the 
functions (pu and 0g. For instance, if one sets (pu — ^ 
then the source fails to excite diffusive momentum den- 
sity in the medium. The full momentum dependence of 
the source is necessary to predict the spectrum and shape 
of the sound wave excited in the medium. The sensitiv- 
ity of the medium's response to the specific mechanism of 
energy and momentum deposition creates the need for a 
hydrodynamic source term derived from first principles. 

In order to derive such a source term one must first 
understand the mechanism of energy and momentum de- 
position and the different scales involved. Consider a 
weakly coupled QCD plasma at an asymptotically high 
temperature T. In this system a fast parton interacts 
with the medium and creates a disturbance at a dis- 
tance scale of the order (gT)^^, where g is the running 
coupling. The disturbance in turn interacts with the 
surrounding medium, creating a new disturbance which 
again interacts with the medium, and the process contin- 
ues, each time at larger and larger distance scales. One 
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could, in principle, determine the medium's response to 
a fast parton at arbitrarily large distances by calculating 
each subsequent disturbance up to the desired scale. In 
practice, however, it is desirable to use an effective the- 
ory, which becomes valid at some length scale, to model 
the medium's response at long distances. As mentioned 
previously, it is believed that nearly ideal hydrodynamics 
is an accurate model of the long distance behavior of the 
QGP produced at RHIC. Hydrodynamics becomes valid 
at length scales which are large compared to the mean 
free path. Again, considering a weakly coupled QCD 
plasma at high temperature, T, the transport mean free 
path is of order {g^T)~^ [2^. Thus to derive an effec- 
tive QCD hydrodynamic source term one should calcu- 
late the medium's response up to distances (at least) of 
order {g'^T)~^, at which point the system evolves hydro- 
dynamically. 

Recently, Neufeld et al. [11] gave a schematic deriva- 
tion of the hydrodynamic source term expected from a 
fast parton moving through a QGP in the relativistic 
limit (velocity approaching the speed of light) in the ab- 
sence of in-medium screening by including the medium 
response at a distance scale of the order {gT)^^ . The 
authors found a propagating Mach cone in the resulting 
linearized hydrodynamic evolution. In this work, I will 
explicitly derive this distribution in the context of kinetic 
theory including color screening. The approach will be 
to treat the fast parton as the source of an external color 
field interacting with a perturbative QGP through the 
coUisionless Boltzmann equation. After obtaining the 
microscopic equation of motion, I find the macroscopic 
evolution by taking moments via the usual Chapman- 
Enskog approach and making the assumption of local 
thermal equilibrium. I find that the fast parton serves 
as a source distribution of energy and momentum for the 
hydrodynamic equations of the medium. I then perform 
an explicit evaluation of the distribution both with and 
without including the medium's response to screen the 
fast parton's color field and interpret the results. 

As mentioned above, a consistent derivation of the hy- 
drodynamic source term should include medium response 
to distances (at least) of order (g'*T)~^. In that sense the 
distribution derived here can be thought of as a sophisti- 
cated representation of a delta function. The calculation 
of the source at intermediate scales will be presented in 
a future publication (in what follows, I choose units such 
that c = ki, = 1). 



bution is given by 

d p „ dp „ dQ" d \ ^, 

a^ + |-^^^ + i-^^+^aQ^)^(-'P'^'^) = o 

(2) 

where the collision term is omitted. The chromodynamic 
equations of motion for a parton with charge gQ°'{t) in a 
color field A'^^ are given by [1^ 



dp 

Itt 



gQ^it) (E'^ + (vxB)") 



and 
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where — dXf^/dt, fabc are the SU{i) structure num- 
bers, and 



+ (v X B)" = 



(5) 



is the color Lorentz force. Inserting ^ and @ into ^ 

yields 



x/(x,p,t, Q) =0. 



(6) 

The singlet and octect components of the parton distri- 
bution, /(x, p,<) and /"(x, p, i), are obtained by taking 
the moments of /(x, p,t, Q) in the color sector: 



/(x,p,t) = j dQ/(x,p,t,0) = /o 



(7) 



and 



P, = j dQ QV(x, p, t, Q) = f1. (8) 

The notation /f used in ([5]) is meant to emphasize that 
any contribution from the color octet distribution must 
come from off equilibrium effects since it vanishes in equi- 
librium. Assuming gA'^^ can be treated as a perturba- 
tion in an otherwise equilibrated system the contribution 
from fi should be small and higher moments, such as f!^'^ 
{a b), which are of higher order in (7^°, will be ignored. 
Applying the first two moments to ^ gives 



(9) 



II. HYDRODYNAMIC SOURCE TERM 



and 



9C2 



(10) 



Consider a system of partons in the presence of an ex- 
ternal color field A° and described by the distribution 
/(x, p, t, Q) which includes the color sector through its 
Q dependence. The Boltzmann equation for this distri- 



where D^ff = d^ff + gfabcA^^fi is the covariant deriva- 
tive. In obtaining (fTO|) I have used 



j dQQ-Q'f{^,p,t,Q) 
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where C2 is the eigenvalue of the Casimir operator of 
the medium partons and, as mentioned above, is ne- 
glected. 

An equation for /q can be obtained by solving pop 
for /f and inserting the result into ([5]). Neglecting any 
space-time dependence in /q the solution is given by 



(12) 



V ■ k + ie 



where 



Uab{x,x') = Pexp 



(13) 



is the path ordered gauge connection. The result 
can be simplified with a contour integration in . It has 
one pole which is in the lower complex plane (i.e., t' < t). 
I find 



/r 



gC2 



d x'Uab{x,x') 
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7V2 - 1 J (27r)3 
^gw.k(t'-t)-,:k.(x'-x)pb(^/^ . ^^j^ 

- f dt'Uab{x,x')F\x')-Vpfo. 



(14) 



where now x' = x(t') = v{t' — t) + x. Combining (fT4|) 
with ^ gives the result for /o [2y| 



(15) 



where 



D, 



N? - 1 



dt'Ftix)Uabix,x')Fhx'). (16) 



The result pS]) describes the phase space distribution 
of partons in the presence of a soft external color field A° 
(the term soft implies the momenta in A^^^ are small com- 
pared to the average momentum in /o). My approach 
is to consider the external color field to be generated by 
a fast parton propagating through the medium, which I 
specify to be a perturbative QGP. In light of the pertur- 
bative nature of the medium it makes sense to expand 
the path ordered gauge connection in p6p to zeroth or- 
der in g. In the hard thermal loop (HTL) approximation 
the source field for a fast parton with constant velocity 
u at position x = ui has the Fourier representation 



EJ(a;) 



(27r)3 

igQp 

(27r)3 



^4^ (ufc' ~ k(k ■ u)) ujSjuj - u ■ k) 



F(fc2-ej,(k,w)tj2) 



^4^g-»fc.x(k-u)(5(w-u-k) 



ei(k, oj)Ljk'^ 



(17) 
k (18) 



B"(a;) 



igQp 
(27r)3 



d'^k 



^-ik-x i^^ u)(5(u;-u-k) 
k"^ — eT(k, Lli)uj'^ 



(19) 



where fc^ ~ (fj^jk), and the dielectric functions, ei(k, w) 
and eT(k, uj), encode the medium's response to screen the 
fields. These functions read explicitly [l^ 
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eT = 1 
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k — UJ 

2\ 



(20) 



where 



/a/3 and Q{k^ — w^) is a step function 
which is unity if k^ > uj'^ and zero otherwise. The deriva- 
tion and momentum space representation of Eqs. (|17l - 
[Tg)) is given in Appendix A. Note that in ((T71) and p5|) 
the electric field has been divided into transverse and 
longitudinal parts such that E — 4- E^ . 

The macroscopic equations of motion arc found by tak- 
ing moments of the evolution equation (|15p : 







(21) 



where = {E, p). The integrals in (PT|) can be evaluated 
by boosting to the frame comoving with the volume ele- 
ment and then exploiting the hydrodynamic assumption 
of local thermal equilibrium by using the thermodynamic 
relations (which are only valid in the local frame of the 
medium element) 



dp p^p" 
(27r)3^ 



/o = (e+p)5"^<5' 



pg 



fiiy 



(22) 



where p is the local pressure, e is the local energy density 
and g'"' is the metric tensor. Using the notation 



(27r)3 ^ ^PiJ) - 



(23) 



I find that the resulting equations of motion for the 
medium evolution are 



where the energy-momentum tensor is given by 
T'^'' = (e + p)w^u''' - pg^""" 



(24) 



(25) 



and is the medium 4- velocity. 

The result ([M]) describes the macroscopic evolution of 
a perturbative QGP in the presence of J'^ , which acts 
as a source of energy and momentum. In other words, 
J" gives the distribution of energy and momentum de- 
posited into a perturbative QGP due to the passage of 
a fast parton. This is precisely the quantity I set out to 
evaluate. In the remainder of this work, I will evaluate 
the source term both with and without dielectric 
screening and examine the resulting distribution. 
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III. EXPLICIT EVALUATION OF THE SOURCE 
TERM 



I begin by simplifying (|23p with an integration by parts 

(26) 

Po 



(2^) 



where 5qi, is the Kronecker-Delta symbol and, as men- 
tioned previously, j range from 1 to 3. I specify the 
medium to be a locally thermal plasma of masslcss gluons 
with the distribution 



/o 



2{Nl 



1 



S/3p" - 1 



(27) 



where Nc is the number of colors and = T is the local 
temperature of the medium. With this specification the 
distribution, /q, now contains the only dependence upon 
the magnitude of the momentum, p'^ ~ p, in (|26p . It 
follows that 



(27r 



Po 



47r 



((^oi/Vj + (5i,.)-D' 



(28) 



where v = p/p" is the directional velocity of medium 
particles. 

Inserting the explicit form of Dij gives 



r{x) 



^ — ik-x—ik' 



dt' I d^kd^k'dw^x 
An 



(29) 



where ttid = \/ Nc/3gT is equal to the Debye mass for 
gluons in the HTL approximation, and I have introduced 
the Fourier representation of the color Lorentz force 



d'^ke-'^-''F^{k). 



(30) 



The generalization to include Nf flavors of 
quarks/antiquarks in the medium can be done in a 
straightforward way by adding the appropriate thermal 
distribution to (P7)) 



2(iV2 - 1) ^ m^Nf 



(31) 



and repeating the steps which lead to ([29|) . The 
only modification is that the expression for the Debye 
mass appearing in (|29p now includes flavor: mo = 
\J {2Nc + Nf)/6 gT. In the numerical results that follow 
I have set Nf = 0. 

I next make the assumption that the typical momen- 
tum transfer between the source particle and the medium 
is small compared to the typical momentum of a medium 
particle. This assumption allows a straight line approxi- 
mation relating the position of a medium particle at time 
t' to its position at time t: 



(32) 



x' =x + v(t'-t). 



The t' integration can be evaluated after plugging p2| 
into (|29p and introducing an exponential damping factor 
ie, yielding 



(27r)8 



V ■ E°(fcO ((5o.v ■ E°(fc) + S,, {Efjk) + (v X B)f (fc))) 
47r(w' - k' • V -f k) 



(33) 



Next consider the integral over v. It can be evaluated 
by choosing a frame in which k' — k'z, performing the 
integral, and then rotating back into a frame in which 
k' is arbitrary. The explicit details of the calculation are 
given in Appendix B. The result is most easily expressed 
by dividing the source term into two pieces, and J2, 
where J'^ + J2 = J'^ ■ After integrating out the delta 
functions in ([l7]-[Tll) I find 

Jr(x) =^ J d3fcd3fc'eHk+k') (-uO^^(k')(^o.E'^(k') 
•E''(k)+5,,(E'^(k') X B''(k)),) 

(34) 



and 



(27r)» J 

X (E''(k') • k') (^(r2i(k')u • k' - (k' 
(5o.(E'^(k).k')+(5..(k' xB'^(k)), 



where 



k' + u-k' 



u k' 



2k''- 



(35) 



iTTU ■ k' 1 

—J^ - T7 (36) 
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In 



k' + u- k' 



{k" - (u • k') 



l\2\ 



k' - U - k' 

- (u • k')^) + 2fc'u • k') 



(37) 



In the above expressions k' is the unit vector in the k' 
direction, k' is the magnitude of k', u is the source par- 
ton's velocity. The exponential damping factor in (|33p . 
le, sets the sign of the imaginary terms in (|36p . (j37p . 

At this point it is convenient to specify a direction 
for the source particle's velocity which I do by choosing 
u = M z. I also choose to work in plane polar coordinates, 
kr and 0, such that kx = kx cos (j) and ky = fc-r sin (/). 
With these choices the only dependence upon the mo- 
mentum space variable (j) is in the exponential and terms 
proportional to cos (j) or sin (j). The exponential depends 
upon (j) through the term ikx (x cos (f> + y sin 0) which can 
be rewritten as ikTpcos[(j)— a], where x = pcosa and 
y = psina. With these simplifications the entire (j) de- 
pendence of and can be integrated out by using 
the relations 



2tt 



2tt 



exp [ikxpicos [(f> — a])] 



1 

cos ( 
sin f 



Joipkr) ' 
fJiipkr) 



(38) 



where J,i(x) is a Besscl function of order n. The cal- 
culation is further simplified due to the fact that the k 
and k' dependence can be completely factorized, reducing 
the remaining four-dimensional integration to a product 
of two-dimensional integrations. Using the explicit field 
dependence listed above, I find there arc 12 unique two- 
dimensional integrations which must be performed. They 
are given by 



.^12. 



dkzdkr 



— I dkzdkx 



kDL{k) 



JoipkT)kTk^' 



dkzdkr- 



— / dkzdkr 



k^Drik) 

^i{z—ut)kz 

DL{k) 

^i{z—ut)k:, 

k^Drik) 



jQ{pkT)kTkl 

Jo{pkT)k^ 
Ji{pkT)kzkr 

■h{pkT)kzk^ 

■T)kl 

,2 

Jo{pkT)kzkT 

'Ji{pkT)k^k'^ 
Ji{pkT)k^ 



Ji{pkrr)k^k^ 
Ji{pkT)krp 



Ji{pkT)k'^kl 
Jo{pkT)kzk^, 



(39) 



where ^^(fc) ee k^eL{k), Drik) = k^ - u'^klerik). 
In terms of the above definitions the source term is 



g^T^iQi? x{y) 
167r4 p 

X (w^Cseg + (^1^7 



(40) 



'aai)) 



and 



167r4 



167r4 

u«6+e4)(e7 + w2ai)) 



(41) 



(42) 



The problem is now reduced to the evaluation of the 12 
a terms. Before considering the full results of such an 
evaluation it is useful to calculate the above expressions 
without including the medium's response to screen the 
fields. I perform this calculation in the next section and 
obtain an analytical result. 



IV. THE UNSCREENED SOURCE TERM 

The unscreened source term is found by setting the 
dielectric function given in (PP]) to unity, or, equivalently, 
setting DL{k) = fc^ and Drik) = k^ - u^kl in It is 
easy to verify that with this simplification and 
can be rewritten as 



J 



^^^^-^^^i|^^6(«^a + (a-«^6)), (43) 



-«2-^^'e5)(e8-^^'a2)) 



A quick perusal of the source term given by (|41l43l44p 
reveals that the only combinations of •■•j ^2 that need 
to be evaluated are (^ — ^^^12) and <^g. These relevant 
combinations are given by 



dkzdkx 



^i{z — ut)k~ 

kl + 72^2 



■y^Ji{pkT)k^ 




^9 


Jo{pkT)kzkT 




a - u^^i2 



(45) 

where 7^ = (1 — u^) ^. 

The above expressions can be evaluated in a straight- 
forward manner by first performing a contour integration 
in the k^ variable. The general form needed is 



dk. 



^±ik~ (z—ut) 



{kl + /fc2 



'I'kl 
kz 



ire 



^{z—ut)kT'y 



1 

±i 



(46) 



where ± refers to the sign of (z — ut). Similarly, one can 
use the relation 



dkrkre-^-^^^'-'''^ 
1 



Jiipkx) 
Joipkr) 



(/?2+72(z-ut)2)3/2 

to obtain the result 

7r7 



P 

7I2 — ut\ 



(47) 



a 

$8 - u'a 



(p2+^2(^_yi)2)3/2 



P 

i(z — ut) 



(48) 
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Before attempting an explicit evaluation of , . . . , it 
is worthwhile to consider the general form of the source 
term given by (|41|43l44p . It is interesting to note that 
(Plj) and ([33]) have a similar form. If one were to re- 
place ^1 — > (u^^2 — w^^s) in ((4T|) then one would find the 
relation 



J, -> uJq 



(49) 



I am motivated to find a relationship between , ^2 , and 
^5 that would allow me to write the source term in a form 
similar to that given by I can (trivially) write 

a = (^"6 - u^a) + (a - ^'6) + (50) 

which allows 

J = 4 [u a <6 + ii) - [[u a - u a) 

(51) 



and 



riv) ^ 2x(j/)^^^7 ("'a + (6 - u%)) (52) 



where the function d is given by 



d{p,z,t) 



87r(p2 4-^2(^_yi)2)3/2- 



(53) 



Notice that J^^J*) now has the same form as — uJ'^ 
apart from the factors of x and {z — ut). The source 
term can be written in the following compact form: 



where 



r{x) = (J"(x),uj°(x) - Jv) 



27 



(54) 



Jv = (x - ut) ^d{p, z, t) {u^^5 + (a - "'6)) • (55) 

I now continue by examining the terms ^7 ■■■iS.6 which 
are considerably more difhcult to evaluate than •■•j C12 
because of their dependence on the two Qi terms. The 
result and method of their evaluation is discussed in de- 
tail in Appendix C. Inserting the ^ into ([331) and ([55)1 I 
eventually find for the source term 



J°(p, z, t) =d{p, z, t)'-fu^ 

X I 1 - ^^7- ^ I 2_ + 



"fup 



{zl 4- p2) 



p2 + 2;2 ^2 



3^{p,z,t) =(x- ut) 



Sir {zl + p2) 

( U^p^ + (z2 ^2 ^ I 2z2 + 



,2^ (9.2 ^ ("'+2)P' 
^2 



V 



(^272 + ^2)^ 



2uz- 



(57) 

where z_ = (z — ut). 

Equations ([55)) and ([57)) . combined with ([M)) . give the 
distribution of energy and momentum deposited into a 
pcrturbativc QGP per unit time by a parton with con- 
stant velocity u = w z at position x — utz, in the absence 
of screening. That this result can be expressed in rela- 
tively simple closed form expressions is remarkable. The 
result for J°{p,z,t) is plotted in Fig. [Tjfor two different 
values of 7. The total energy deposited per unit time 
(or length) can be obtained by integrating ([56)) over all 
space: 



dE 
dx 



dxj"{p,z,t) 



In 



72^ 



(58) 



where y = cosh" "'^(7) is the rapidity and pmax and pmin 
are infrared and ultraviolet cut-offs for the p-integration. 
This result matches the standard leading-logarithmic 
result [13 for collisional energy loss in the relativis- 
tic (u —f 1) limit if one chooses pmax = I/'tt-d and 
Pmin = l/(2^£'pT), where Ep is the energy of the fast 
parton. It is interesting to plot the u dependence of ([55)) 
which is done in Fig. [2] The total momentum deposited 
per unit time can likewise be obtained by integrating 
3{p, z, t) over all space (recall that J = uJ° — Jv)- One 
can verify explicitly that the source term satisfies the 
energy-momentum relation 



J dxu ■ 3{p, z,t) ~ J dxj'^{p, z,t). 



(59) 



A careful observer may notice that the denominators of 
([56)1 and ([57)) contain terms involving both [z^'f^ + p^), 
which encodes the Lorentz contraction of the fast par- 
ton's field configuration, and i^z"^ + p^) , which is sensi- 
tive to the rest frame of the medium. One may ask why 
the result, in the absence of dielectric screening, should 
depend on the rest frame of the medium. The answer 
can be traced back to the t' integration in ([^^ which 
records the history of the fast parton's interactions with 
the medium. In that expression the fast parton inter- 
acts with a medium particle at t' , which is expressed by 
Fj{k'), and again at t, which is expressed by Fi{k). In 
between t' and t the medium particle propagates with a 



7 



for 7 = 20 




-2 




FIG. 1: (color online) Three-dimensional plots of the scaled energy density deposited per unit time by a fast parton moving 
in the positive z direction for Os = I/tt and Nc — "i [see ((55)) ]. The plots are scaled by mnCQp)^ « 20.4 GeV/fm* for a gluon 
moving in a QGP at 200 MeV. As one might expect, the distribution becomes Lorentz contracted for large values of 7. (In the 
above plots, I have set y — 0.4 GeV^^ to avoid plotting the origin where the source is divergent.) 



1- 



cosh ^(7) 



Atu^) with p = 1 [GeV"^] 



u y 




FIG. 2; Plot of the velocity dependence of —dE/dx as given 
by (|58[) where it is in units of c = 1. 



momentum that depends on the properties of the medium 
[see ([32]) ]. Thus the second interaction, at i, depends 
on both the field configuration of the fast parton and 
the properties of the medium, even in the absence of 
screening. One can observe this effect in Fig. |3l where 
J°(p, z, t)/{'-fu^) is plotted as a function of 72 for fixed p. 
If J°{p,z,t) depended only on the field configuration of 
the fast parton then the curves shown in Fig. [3] should 
be identical for any velocity. However, in the case of a 
highly relativistic source parton one sees a slight shift 
in the curve reflecting the dynamics of the medium in 
between the two interactions at t' and t. 



For a highly relativistic source parton Eqs. ([55)1 and 
57)) can be written in a simplified form by considering 



that terms with (2^7^ + p^) in the denominator will be 
severely damped except for z w ut. In this case, I keep 
(z — ut) to first order (and drop terms proportional to 




Y (z-ut) [GeV ] 



FIG. 3: (color online) A plot of [p, z,t) / {-yv?) as a function 
of 7(2 — ut) for the parameters u = 0.99, ttid = 0.4 GeV, 
Ua = I/tt and p = 1 GeV~^. If J°{p,z,t) depended only on 
the field configuration of the fast parton the curves should be 
identical. As one can see, there is a slight shift in the curve 
when going from 7 = 1.25 to 7 = 20. This shift reflects the 
dynamics of the medium in between the two interactions at 
t' and t (see discussion in the text). The fact that the shift 
in the curve is small even though 7 has changed by a factor 
of 16 suggests the effect is small. 



1/7, etc.) and the result becomes 



J°{p,z,t) = d{p,z,t)ju^ 1 



3v{p,z,t) = {x. — ut) d{p, z,t) 



7u z_ 



z^7^ + /O^ 



^17^ + p^ 



(60) 
(61) 



In the next section, I consider the numerical evaluation 
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FIG. 4: (color online) Plots of the scaled energy density deposited per unit time by a parton moving in the positive z direction 
with velocity u = 0.99 and Qs = I/tt, both with and without medium screening (compare with Fig. [l]). There is a very similar 
structure in both cases; however, as one would expect, the screened result shows some damping in the transverse direction. 




FIG. 5: (color online) Plots of the scaled transverse momentum density deposited per unit time both with and without medium 
screening included for the same parameters as in Fig. 2] The structure shows noticeable differences in the two cases. In 
particular, the screened result has a positive bump just behind the source parton which is absent in the unscreened result. 



p = 1.0 [Ge\rM 

Screened 

Unscreened — 



[GeV^] 







-3 -2 -1^ 

-o.ooa 




2 


3 


-0.004| 








-0.006 









(z - u t) 



FIG. 6: (color online) Plot of transverse momentum density 
deposited per unit time by a gluon both with and without 
medium screening included (compare with Fig. [S}. The plot 
is for It = 0.99, ttid = 0.4 GeV, and Qs = I/tt. 



of the source term with medium screening included and 
compare to the analytical result without screening. 

V. RESULT WITH SCREENING 

The evaluation of the source term (|40l41l42p with 
medium screening included must be done numerically. 
As mentioned before this amounts to evaluating the 12 

terms listed in . A discussion of the specific tech- 
niques used to perform the numerical integrations is given 
in Appendix D. 

In this section all evaluations are for a parton mov- 
ing along the positive z axis with velocity u = 0.99 and 
I have set mD = 0.4 GeV, A'c = 3 and = I/tt (the 
choice of w = 0.99 corresponds to 7 « 7, or « 3 GeV 
for the source gluon, which is likely too small a choice 
for the fast partons produced at RHIC but lends itself 
nicely to plotting the result). The result for the energy 
density deposited per unit time both with and without 
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medium screening is given in Fig. [4] One sees a very 
similar structure in both cases except that the screened 
result is moderately damped in the transverse direction. 
A similar plot is shown in Fig. Ofor the transverse mo- 
mentum density deposited per unit time. Here one finds 
more noticeable differences in the two cases. For instance, 
the screened result has a positive bump just behind the 
source parton which is absent in the unscreened result. 
This would correspond to matter being pushed outward, 
away from the source parton, in the vicinity of the bump. 
The differences are further highlighted in Fig. [HI Plots 
of the parallel momentum density distribution, J^, are 
not included because they are very similar to J° for a 
relativistic source. 



VI. DISCUSSION 

In this paper, I have derived the space-time distribu- 
tion of energy and momentum deposited by a fast parton 
traversing a perturbative quark-gluon plasma by includ- 
ing physics at the Debye screening scale. The approach 
has been to treat the fast parton as the source of an exter- 
nal color field perturbing the medium. The final result 
depends on two parameters: the source strength asQp 

and the Debye mass tod = y/Nc/3 gT. 

Several assumptions were made in deriving the result. 
The first, and primary, of these is that the medium is 
perturbative in the coupling as. This assumption is re- 
flected in the derivation of p5|) , the expansion of the path 
ordered gauge connection, Uabix, x'), and the choice of 
fields pn - [T9|) . Another assumption is that the medium 
is in local thermal equilibrium. This is the basic assump- 
tion of hydrodynamics and is necessary to obtain (|24|1 . I 
also assume that the typical momentum transfer between 
the source particle and the medium is small compared to 
the typical momentum of a medium particle. This as- 
sumption allows the straight line approximation given 
by ((32)) . To get an idea of where this last assumption 
breaks down one can use the equipartition theorem for 
masslcss particles to find a relation between the typical 
momentum of a medium particle and the temperature: 
p = ST. For T = 350 MeV the typical momentum of a 
medium particle is about 1 GeV which places an upper 
limit on typical momentum transfers. It follows that the 
straight line approximation breaks down somewhere on 
the order of fc > 1 GeV~^ (or distances less than 0.2 fm). 

In the absence of the medium's response to screen the 
fields, it was possible to obtain an analytical result which 
is given by (j56p and (|57p . This result shows many similar- 
ities to the result obtained when the medium's response 
to screen the color field is included. The most noticeable 
differences appear in the transverse momentum density 
distribution. The total energy deposited per unit time is 
obtained by integrating over all space and is 

given by ([55)1 for the unscreened result. 

The reader interested in reproducing the numerical re- 
sults for the screened source term should consult Ap- 



pendix D, where the specific techniques used to perform 
the numerical integrations are discussed. 
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Energy under Grant DE-FG02-05ER41367. 



Appendix A: Derivation of (|17p - (|19p 

In the HTL approximation, the source fields obey 
Maxwell's equations, which are given by 



V-D(x,<) = p(x,t) 

o»B(x, t) 



V X E(x,t) = -- 
V-B(x,0 = 



dt 



(A-1) 

(A-2) 
(A-3) 



VxB(x,i) = j(x,t) + ^5^Jl^ (A-4) 



dt 



where the source density is 



p(x, i) = gQpS{x-ut) 
j(x,i) = u/9(x,i) 



(A-5) 
(A-6) 



(in the above equations u is the source particle's velocity 
and I have suppressed the color index, which I will restore 
in the final result) . It is convenient to solve for the fields 
by taking the Fourier transform of Maxwell's equations. 
Using the general rule 

= (2^)^ / / ^''"""^"'^e^' ^) 
it is found that 

ik-D(k,cj) = p(k,Lj) (A-8) 

ikxE(k,tj) = iujB(k,uj) (A-9) 

ik-B(k,w) = (A-IG) 
zkxB(k,w) = j(k,u;) - iu;D(k,tj) (A-11) 

where 

p(k,w) = 2TrgQpS{uj-u-k) (A-12) 

j(k,c^) = up{k,Lo). (A-13) 

The function D is related to E through the dielectric 
tensor e,, : 



D, = e,,E, 



(A-14) 



(the summation over j is implied and I now suppress the 
(k, a;) notation). In an isotropic medium the dielectric 
tensor can be decomposed into transverse and longitudi- 
nal parts such that 



fc2 



fc2 



(A-15) 
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and D is now written as 

D = elEl + etEt 



(A-16) 



where El = k(k-E) and Er = E — El. One should note 
that k • Et = and k x El = 0. Maxwell's equations 
now take the form 

zcLk-EL = p (A-17) 

ik X Et = iiuB (A-18) 

ik-B = (A-19) 

ik X B = up - iLu{eLEL + etEt). (A-20) 

An equation for B can be obtained by crossing k into 
(|A-20P which combines with (|A-18p and (|A-19|) to yield 
(with color index restored) 



2nigQ°:(k x u) 
B" = ^S{lu - u • k). 



(A-21) 



One can then obtain Et by crossing k into (jA-lSP 
, 27r*.gO^^ (ufc2 - k(k • u)) 



Combining (|A-21|) and (|A-22|) with (|A-20|) yields 



(A-22) 



i ( k^ 
El = (-1) up + i — Et - iujeT'Er) 



2^igQl{k ■ u) 



(A-23) 



eLwfc^ 



-k(5(tj — u • k). 



These results can be back transformed to position space 
to yield ([17] -[HI). 



Appendix B: Explicit evaluation of the v integration 

In going from Eq. (|33p to and (|55|) I performed an 
integration over v by choosing a frame in which k' = k'z, 
performing the integral, and then rotating back into a 
frame in which k' is arbitrary. In this appendix the ex- 
plicit details involved will be given. The relevant starting 
point is 



V ■ E'" {Sq^v ■ E° + 4. {Ef + (v X B)f)) 
47r(u;' - k' • V + ie) 



(B-1) 



where E'° is short for E''(fc'), etc. There are two basic 
integrals to consider, which have the form 



47r([j' - k' • V + ie) 
v,v„ 



47r(cj' - k' • V + ie) 



(B-2) 
(B-3) 



Now specify that k' = k'z and consider that since 
V = (sin 9 cos cj), sin 9 sin (f), cos 9) the integral over d(j) will 



ensure that only survives in (jB-2p . Making the re- 
placement dv ^ J d9d(l) sin 9 gives 



Airiuj' — k' cos 9 + ie) 
1 ... (B-4) 



da 



(w' - k'^ + ie) 



Similarly, in (|B-3p the integral over dcj) ensures only I^x, 
bute: 

1 d^i-e) 



lyy, and Izz contribute 



and 



1 



1 {uj' - fc'e + ie) 
dU' 



2 y_i (cj' - k'^ + ie) 
The terms are evaluated with the result 

k + LO 







In 


fc -|- a; 


k^' 


k — UJ 


— 7ri(l 


k^' 


+ 2 


1) 



fi3(fc) 



2fc2 V /fc 



In 



k-i 



TTIUJ 
~k' 



Equations (|B-4I - IB-6P can be rewritten as 



-^jrn i^jx^nix ~^ ^jx^mx) ^2 ^jz^mz^S 
= Sjm^2 + SjzSmz{-^^l — ^2)- 



(B-5) 
(B-6) 
(B-7) 

(B-8) 

fliik). 
(B-9) 

(B-10) 
(B-11) 



In the above expressions the subscript z denotes the 
orientation of the k' vector. Generalizing to an arbitrary 
k' is done by making the replacement Sjz — » k'j in (jB-10|) 
and (jB-lip . Going back to (jB-l[) (and suppressing the 
color index) I find 



dv 



V • E'° (<5o.v • E'' + S,, {Ef + (v X B)°)) 



47r(a;' - k' • v -I- ie) 

^OiyEjEjnljrn + {EjEilj + eimkl jmE'jBk) 

So, ((E' • E)Q2 + (k' • E')(k' • E)(^r!i - 
5,. ((k' • E')-B«Oi + 

((E' X B),r!2 + (k' • E')(k' X B),(^r!i - n^) 



(B-12) 



which is the result given in ([M]) and (|35p . 
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Appendix C: Evaluation of ^i, ...,^6 



In Sec. IIVI I put off the explicit evaluation of the terms 
^1, ...,.^6 and directly went to the results and ([57)1 . 
In this appendix, I will give a detailed derivation of the 
expressions: 



= / dkzdkT 



dkzdkT 



e^i^-ut)k,ni{k) 
kDL{k) 

k^Drik) 



'MpkT)kTk'^' 
Jn{pkT)kTkl 

Jo{pkT)k^ 
.Ji{pkT)kzk^. 

Jo{pkT)kzk^ 

,2 1,2 



Ji{pkT)k^k'^ 



(C-1) 



Here, I am interested in the unscreened source term 
which is found by setting DL{k) = k^ and Drik) = 



i^fc^. I start by changing to polar coordinates (i.e., 
k and 9, where kz = kcosO and fcy = fcsin^) in which 
case the two terms can be written as 



Vli{k) 



2k 

G2{e) 

Ak 



(C-2) 
(C-3) 



where 

Gi{0) = 2MCOs6'tanh"^[ucos6'] - inu cos 9 - 2 (C-4) 
and 

G2{9) = 2tanh-^[ucos6l](l-u2cos2 6l) 



— Z7r(l — cos^ 9) + 2ucos( 



(C-5) 



Defining A(0) as some generic function, I now write the 
general form of the k integration as 



= / dkd9 sin 9e 



ik{[z — ut) cos Q) 



Jo{pksui9) 
Ji [pk sin 9) 



which can be evaluated using the relations [28| : 

1 



(C-6) 



dke^'^yjoikb) = 

y^b^ -{y + ie)2 



1 / 

dke'^yji{kh) = 1 



^\P- - (y + ief 




where 6 > and e is a positive infinitesimal quantity. 
Defining a = — ut)/ p and using the above integrals 
allows 



= / dfcd6'sin6ie*'=((^-"*)™''*) 
,1 



Jo{pk sin 9) 
Ji{pksin9) 



A,(0) 



d9- 
P 

d9- 

P 



■\J sin^ 9 — {a cos O^ia)'^ 



1 



p-\J sin^ 9-~{a cos 0-\-ie)'^ 



MO) 



U{9) 

m 



HO)- 



(C-9) 



where I have set t 
terms look like 



for conciseness. Explicitly, the 



d9-G,{9) 
2p 

^^^G^W 



U{9) 
U{9) cos^ 9 
L{9) cos 9 sin 9 

U{9) cos 9 sin^ t 
L{9) cos2 6'sin^ 



(C-10) 



where ^3 is omitted since it will not appear in the final 
expressions. 

It is imperative to consider symmetry before trying to 
evaluate the various terms. For instance, in ^1, ...,^4 the 
imaginary part of U [9) and the real part of L{9) only need 
to multiply against the term —iiru cos 9 due to symmetry 
considerations. Similarly, the real part of U{9) and the 
imaginary part of L(9) do not need to multiply against 
~iTTucos9. It turns out that due to symmetry the terms 
involving the imaginary part of U(9) and the real part of 
L{9) are easier to evaluate than those involving the real 
part of U{9) and the imaginary part of L{9). It makes 
sense to first get these simpler expressions out of the way 
and then focus on the more tedious ones. I introduce 
the notation ■Ci = ?ia + Cifej where now denotes terms 
involving the imaginary part of U{9) and the real part of 
L{9). Using Mathematica 6.0 these terms can be evalu- 
ated and I find 



^5 a 
^6 a 



7rau(3+a ) 

i-rvu / 1 _ z a(3+a 
' p X* 

_ z a(3+a^) 
P X'' 



(C-11) 



6p 



where x = Vl + a^. 

Now I turn to the more involved terms which in- 
volve the real part of U{9) and the imaginary part of 
L{9). It is convenient to note that 



Re[U{9)] = sin 9 Re 



1 



\/ sin^ 9 — cos^ 
1 

a/ sin^ 9 — a? cos^ 1 



(C-12) 
(C-13) 



For the sake of bookkeeping, I note there are 7 distinct 
integrations which need to be performed. These terms 
are given by the real parts of 



Vt= I d9 



smf 



\/ siv? 9 — cos^ 



(C-14) 



where 



























T7. 





1 

0082 6* 

COS 9 tanh^ ^ [u cos 9] 
cos'^ 9 tanh^ [u cos 9] 
cos 9 sin^ 9 tanh~ ^ [u cos 9] 
cos2 6'sin2 6'/ (1-^2 cos2 9) 
cos'' 9/{l-u^ cos^ 61) 



In terms of these the are 









1 








P 


-$66- 





Uip3 - (^1 
U(^4 - (^2 

((y55 + uipe)/2 



(C-15) 



(C-16) 



Using Mathematica 6.0 1 find for the real parts of the ip 
terms: 



(C-17) 



= — J ( X - + 1/7^ 



(X^ + 2^2) 7i_ 



7%4 \ 



TT / 1 1 
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(C-19) 

(C-20) 
(C-21) 

(C-22) 



(C-23) 



tP2 



2x3 



(C-18) Combining the above results I find 
I 



^5 



-TT ^-^/a^ + 1/7^ — au^ 



TT (2^a^ + ^^72 (a3«3 + x^) + (1 + a^^^) (2u2 + x2(l ~ Su^)) - Sx^) 



(2 + + x^ (x - 3a^^) + 2ya2 + ^(x^ - u')^ 



u6 - W ^5 



u^e + 2^4 



TTM ^^a (^-207^ + 2^72 y'a^ + + a) - 1 



2x%/a2 + ^72 



(^-a2 + (a^ + 3) ua^a"^ + + 2u'^ - 1^ 
4 /„2 I 1 



ITTU 



(C-24) 



which are combined with ([48|) to yield ([56|) and ([57]) . grations. The terms are given by 



Appendix D: Numerically Integrating Equation (|39p 

As mentioned in Sec|Vl the evaluation of the 12 
terms listed in (|39p with medium screening included must 
be done numerically. In this appendix I will discuss the 
specific techniques used to perform these numerical inte- 



— \ dkzdkT 



kDL{k) 



Jo{pkT)kTkl 

JoipkT)k^ 
.Ji{pkT)kzkj 



(D-1) 
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dkzdkx 



JoipkT)kzk^ 
Ji{pkT)k?^kl 



— I dkzdkx- 



DL{k) 



dkzdkx 



^i{z—ut)k- 

k^Drik) 



Ji{pkT)kl ' 
Jo{pkT)kzkT 



Ji{pkT)k'^k'^' 

Ji{pkT)k^ 
Ji{pkT)k'^k'l 
.JoipkT)kzk^, 



(D-2) 
(D-3) 

(D-4) 



where the k^ integration runs from ±00 and the kx in- 
tegration ranges fr'om — > 00. Since each S,i contains a 
double integration which ranges over an infinite interval, 
care must be taken to cast each term in a form that will 
tend to zero as quickly as possible in order to make nu- 
merical evaluation feasible. In the case of the integrals 
over dkz , we can take advantage by the exponential term 
gi(z-iit)fej |-jy Jjgriding the contour in the imaginary plane. 
The strategy works as follows: start with a function of 
the form 



-^*^'''f{kz)dkz 



(D-5) 



where f{kz) is an analytic function. Next, bend the con- 
tour of the integral in the complex plane by making the 
variable change 



and 



dk, dre'"S'sn[r] 



(D-6) 



(D-7) 



where the limits of integration on r are from —00 to 00. 
Under this change, (|D-5P becomes 

cxp [i{z - wt)re''"S's"M]/(re'"S's"M)e*"'^'snW^^ 

_ J 

(D-8) 

and now the exponential term is damped in (|D-8p [pro- 
vided a has the same sign as (z^ut)] making the integral 
more numerically manageable. (|D-8[) and (|D-5p will be 
equivalent if the integrand vanishes as r — > ±00 and no 
poles from f{kz) are crossed when bending the contour 
(this places a restriction on the value of a). 

Before applying this technique it is necessary to con- 
sider where the have poles. Any poles located directly 
on the imaginary axis will not affect the result so we 
can ignore those. That leaves us to consider where the 
zeros of Dxiukz) and DL{ukz) are. An analytic solu- 
tion for the zeros of Dxiukz) and D^iukz) is impossible 
because of the logarithmic terms. However, the zeros 
can be searched for numerically. Using the parameters 
u — 0.99z and ttt-d = 0.4 GcV, it is found that all poles 
are sufficiently close to the imaginary axis so that the 
value of a used is not restricted for any reasonable value 



for terms with Driukz) in the denominator. I choose to 
use |a| = 1. However, for terms with DL{ukz) in the de- 
nominator, it is found that for values of fc^ ^ 1/2 poles 
become increasingly close to the real axis. For this rea- 
son, I choose to break up the integration in the terms 
with DL{ukz) in the denominator into a term 



dkx 



dk. 



(D-9) 



which I bend in the complex plane with a value of |a| = 1 
and a term 

dkr / dkz (D-10) 

J —oc 

which I do not bend in the complex plane. For the terms 
I do not bend in the complex plane it is usually necessary 
to subtract off the asymptotic form of the kz integration 
in order to make the numerical integration efficient. The 
asymptotic form must then be added back. Similarly, for 
the terms I do bend in the complex plane, I usually must 
subtract off the asymptotic form of the fc^ integration to 
provide efficient convergence. 

In the case of z = ut, it does not make sense to bend 
the kz integral in the complex plane since there will be 
no damping term from the exponential. In this case, I 
instead perform the numerical integration by subtracting 
off the unscreened form of in the integrand (which has 
the same large k form) and then re-adding after evalua- 
tion. 

To make these ideas more concrete, consider ^7 for the 
case oi z ^ ut. I write 



dk. 



dkre'^''-'"'^^' Jiipkr) 



1 



DL{ukz 

(D-11) 

where I have subtracted 1 to make the kx integration 
more efficient. I bend the above contour in the complex 
plane and perform the integration and then re-add 

dkz I dfcTe'(^-"*)'=-Ji(pfcT)-(5(z-Mi) (D-12) 
-00 Ji 

which can be ignored, since it only contributes when z = 
ut. Next, I integrate 

dkz C dkre'^'-^'^'^' Ji{pkT) ^ f, , 
-1 Jo DL[ukz) 

1 / 7,2 



(D-13) 

where I have subtracted k^ / k1 to make the kz integration 
more efficient. I then must re-add 
»i 



2Re / dkz / dkre'^^-""'^^^ .h{pkT)-^ 

Jl Jo . 

Up) 



(D-14) 



(2(cos(z)-t-zSi(z))-7r|z|). 
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Considering ^7 for the case oi z = ut, I write 



and evaluate at z == ut. 



/oc poo 
-00 Jo 



(D-15) 

and then analytically re-add [see (|48)) ] 



'^T + «3 (D-16) 







(p2 + (z _ Mt)2)3/2 
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Fast Partons as a Source of Energy and Momentum in a Perturbative Quark-Gluon 
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I derive the space-time distribution of energy and momentum deposited by a fast parton traversing 
a perturbative quark-gluon plasma by considering the fast parton as the source of an external color 
field interacting with the medium. I include the medium's response to screen the fast parton's color 
field by incorporating dielectric functions and compare to the unscreened result. 



PACS numbers: 12.38.Mh 

I. INTRODUCTION 

Quantum chromodynamics (QCD) predicts a transi- 
tion in sufRciently hot and/or dense nuclear matter from 
colorless hadrons to a state of deconfined quarks and glu- 
ons known as the quark-gluon plasma (QGP) [H- It is 
believed that this transition has been observed in ex- 
periments performed at the Relativistic Heavy-Ion Col- 
lider (RHIC) (see, for example, [1]). Fast partons (par- 
tons with velocity approaching the speed of light) created 
by hard transverse scattering in the early moments of a 
heavy-ion collision are a useful probe in understanding 
QGP physics. Of particular importance is the process 
of jet quenching in which fast partons lose energy and 
momentum by interacting with the surrounding medium 
(see, for instance, I, i, i, S 0, S i, ) • 

While much effort has been devoted to understanding 
the dynamics of jet quenching, the question of how the 
energy and momentum deposited by the fast parton af- 
fects the bulk behavior of an evolvin g Q GP is relatively 
new (see, for example, [11 [11 [11, QIll [11, [13, EH)- 
One promising approach is to treat the fast parton as the 
source of energy and momentum coupled to the hydrody- 
namic equations of the medium. There is strong evidence 
[20I . [2l| that the matter created in heavy-ion collisions at 
RHIC can be described by nearly ideal hydrodynamics 
and it was argued in [Tsj that if the energy and momen- 
tum deposited by a fast parton thermalize quickly the 
resulting disturbance should evolve hydrodynamically as 
well. 

Casalderrey-Solana et al. [ij] suggested that since the 
fast parton travels through the medium at supersonic 
speed, one might expect a Mach cone to develop in the 
surrounding hydrodynamic medium. They proceeded to 
show that if one uses an appropriately chosen supersonic 
sound source in a linearized hydrodynamic evolution one 
indeed observes a propagating Mach cone. In their work, 
the explicit calculation of the source term generated by 
a supersonic parton was not addressed. In a later work 
I13I Chaudhuri and Heinz examined a two-dimensional 
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hydrodynamics simulation of an expanding QGP with a 
phenomenologically motivated ansatz as their fast parton 
source term. Here the authors failed to find the formation 
of a Mach cone except for seemingly unphysical values of 
energy loss. Much recent work [ll, [ll| has been focused 
on examining a supersonic color charge passing through 
an A/" = 4 supersymmetric Yang-Mills plasma (in both 
of the above citations the authors find a Mach cone in 
the resulting dynamics). Such efforts are motivated by 
the idea that the system is a useful model of a real QCD 
plasma. 

As the above examples illustrate, there have been sev- 
eral insightful attempts at understanding the influence of 
fast partons on the bulk evolution of an expanding QGP. 
However, the explicit evaluation of the distribution of 
energy and momentum deposited by a supersonic parton 
traversing a QCD plasma has remained an open ques- 
tion. It has been pointed out [ll|, [l^l that the medium's 
response to a fast parton is sensitive to the specific form 
of energy and momentum deposition. For example, the 
hydrodynamic source vector of a parton moving at con- 
stant velocity, u, can be written in Fourier space as (Toj 

J(cj,q)= [vL(j)u{uj, q^) +iq(t)q{uj,q^)\2TT6{uj -u-ci). (1) 

The medium's response then becomes sensitive to the 
functions (pu and 0g. For instance, if one sets (pu — ^ 
then the source fails to excite diffusive momentum den- 
sity in the medium. The full momentum dependence of 
the source is necessary to predict the spectrum and shape 
of the sound wave excited in the medium. The sensitiv- 
ity of the medium's response to the specific mechanism of 
energy and momentum deposition creates the need for a 
hydrodynamic source term derived from first principles. 

In order to derive such a source term one must first 
understand the mechanism of energy and momentum de- 
position and the different scales involved. Consider a 
weakly coupled QCD plasma at an asymptotically high 
temperature T. In this system a fast parton interacts 
with the medium and creates a disturbance at a dis- 
tance scale of the order (gT)^^, where g is the running 
coupling. The disturbance in turn interacts with the 
surrounding medium, creating a new disturbance which 
again interacts with the medium, and the process contin- 
ues, each time at larger and larger distance scales. One 
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could, in principle, determine the medium's response to 
a fast parton at arbitrarily large distances by calculating 
each subsequent disturbance up to the desired scale. In 
practice, however, it is desirable to use an effective the- 
ory, which becomes valid at some length scale, to model 
the medium's response at long distances. As mentioned 
previously, it is believed that nearly ideal hydrodynamics 
is an accurate model of the long distance behavior of the 
QGP produced at RHIC. Hydrodynamics becomes valid 
at length scales which are large compared to the mean 
free path. Again, considering a weakly coupled QCD 
plasma at high temperature, T, the transport mean free 
path is of order {g^T)~^ [2^. Thus to derive an effec- 
tive QCD hydrodynamic source term one should calcu- 
late the medium's response up to distances (at least) of 
order {g'^T)~^, at which point the system evolves hydro- 
dynamically. 

Recently, Neufeld et al. [11] gave a schematic deriva- 
tion of the hydrodynamic source term expected from a 
fast parton moving through a QGP in the relativistic 
limit (velocity approaching the speed of light) in the ab- 
sence of in-medium screening by including the medium 
response at a distance scale of the order {gT)^^ . The 
authors found a propagating Mach cone in the resulting 
linearized hydrodynamic evolution. In this work, I will 
explicitly derive this distribution in the context of kinetic 
theory including color screening. The approach will be 
to treat the fast parton as the source of an external color 
field interacting with a perturbative QGP through the 
coUisionless Boltzmann equation. After obtaining the 
microscopic equation of motion, I find the macroscopic 
evolution by taking moments via the usual Chapman- 
Enskog approach and making the assumption of local 
thermal equilibrium. I find that the fast parton serves 
as a source distribution of energy and momentum for the 
hydrodynamic equations of the medium. I then perform 
an explicit evaluation of the distribution both with and 
without including the medium's response to screen the 
fast parton's color field and interpret the results. 

As mentioned above, a consistent derivation of the hy- 
drodynamic source term should include medium response 
to distances (at least) of order (g'*T)~^. In that sense the 
distribution derived here can be thought of as a sophisti- 
cated representation of a delta function. The calculation 
of the source at intermediate scales will be presented in 
a future publication (in what follows, I choose units such 
that c = ki, = 1). 



bution is given by 

d p „ dp „ dQ" d \ ^, 

a^ + |-^^^ + i-^^+^aQ^)^(-'P'^'^) = o 

(2) 

where the collision term is omitted. The chromodynamic 
equations of motion for a parton with charge gQ°'{t) in a 
color field A'^^ are given by [1^ 



dp 

Itt 



gQ^it) (E'^ + (vxB)") 



and 



(3) 



(4) 



where — dXf^/dt, fabc are the SU{i) structure num- 
bers, and 



+ (v X B)" = 



(5) 



is the color Lorentz force. Inserting ^ and @ into ^ 

yields 



x/(x,p,t, Q) =0. 



(6) 

The singlet and octect components of the parton distri- 
bution, /(x, p,<) and /"(x, p, i), are obtained by taking 
the moments of /(x, p,t, Q) in the color sector: 



/(x,p,t) = j dQ/(x,p,t,0) = /o 



(7) 



and 



P, = j dQ QV(x, p, t, Q) = f1. (8) 

The notation /f used in ([5]) is meant to emphasize that 
any contribution from the color octet distribution must 
come from off equilibrium effects since it vanishes in equi- 
librium. Assuming gA'^^ can be treated as a perturba- 
tion in an otherwise equilibrated system the contribution 
from fi should be small and higher moments, such as f!^'^ 
{a b), which are of higher order in (7^°, will be ignored. 
Applying the first two moments to ^ gives 



(9) 



II. HYDRODYNAMIC SOURCE TERM 



and 



9C2 



(10) 



Consider a system of partons in the presence of an ex- 
ternal color field A° and described by the distribution 
/(x, p, t, Q) which includes the color sector through its 
Q dependence. The Boltzmann equation for this distri- 



where D^ff = d^ff + gfabcA^^fi is the covariant deriva- 
tive. In obtaining (fTO|) I have used 



j dQQ-Q'f{^,p,t,Q) 



m - 1 



/o + fr (11) 
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where C2 is the eigenvalue of the Casimir operator of 
the medium partons and, as mentioned above, is ne- 
glected. 

An equation for /q can be obtained by solving pop 
for /f and inserting the result into ([5]). Neglecting any 
space-time dependence in /q the solution is given by 



(12) 



V ■ k + ie 



where 



Uab{x,x') = Pexp 



(13) 



is the path ordered gauge connection. The result 
can be simplified with a contour integration in . It has 
one pole which is in the lower complex plane (i.e., t' < t). 
I find 



/r 



gC2 



d x'Uab{x,x') 



9C2 



7V2 - 1 J (27r)3 
^gw.k(t'-t)-,:k.(x'-x)pb(^/^ . ^^j^ 

- f dt'Uab{x,x')F\x')-Vpfo. 



(14) 



where now x' = x(t') = v{t' — t) + x. Combining (fT4|) 
with ^ gives the result for /o [2y| 



(15) 



where 



D, 



N? - 1 



dt'Ftix)Uabix,x')Fhx'). (16) 



The result pS]) describes the phase space distribution 
of partons in the presence of a soft external color field A° 
(the term soft implies the momenta in A^^^ are small com- 
pared to the average momentum in /o). My approach 
is to consider the external color field to be generated by 
a fast parton propagating through the medium, which I 
specify to be a perturbative QGP. In light of the pertur- 
bative nature of the medium it makes sense to expand 
the path ordered gauge connection in p6p to zeroth or- 
der in g. In the hard thermal loop (HTL) approximation 
the source field for a fast parton with constant velocity 
u at position x = ui has the Fourier representation 



EJ(a;) 



(27r)3 

igQp 

(27r)3 



^4^ (ufc' ~ k(k ■ u)) ujSjuj - u ■ k) 



F(fc2-ej,(k,w)tj2) 



^4^g-»fc.x(k-u)(5(w-u-k) 



ei(k, oj)Ljk'^ 



(17) 
k (18) 



B"(a;) 



igQp 
(27r)3 



d'^k 



^-ik-x i^^ u)(5(u;-u-k) 
k"^ — eT(k, Lli)uj'^ 



(19) 



where fc^ ~ (fj^jk), and the dielectric functions, ei(k, w) 
and eT(k, uj), encode the medium's response to screen the 
fields. These functions read explicitly [l^ 



EL 



1 



eT = 1 



2m^ 

9 

W2 



1 



UJ 

2k 



In 



"a;2 


UJ 


(fc2 


F + 


2k 


1 


X ^In 


'k 










k 





k — UJ 

2\ 



(20) 



where 



/a/3 and Q{k^ — w^) is a step function 
which is unity if k^ > uj'^ and zero otherwise. The deriva- 
tion and momentum space representation of Eqs. (|17l - 
[Tg)) is given in Appendix A. Note that in ((T71) and p5|) 
the electric field has been divided into transverse and 
longitudinal parts such that E — 4- E^ . 

The macroscopic equations of motion arc found by tak- 
ing moments of the evolution equation (|15p : 







(21) 



where = {E, p). The integrals in (PT|) can be evaluated 
by boosting to the frame comoving with the volume ele- 
ment and then exploiting the hydrodynamic assumption 
of local thermal equilibrium by using the thermodynamic 
relations (which are only valid in the local frame of the 
medium element) 



dp p^p" 
(27r)3^ 



/o = (e+p)5"^<5' 



pg 



fiiy 



(22) 



where p is the local pressure, e is the local energy density 
and g'"' is the metric tensor. Using the notation 



(27r)3 ^ ^PiJ) - 



(23) 



I find that the resulting equations of motion for the 
medium evolution are 



where the energy-momentum tensor is given by 
T'^'' = (e + p)w^u''' - pg^""" 



(24) 



(25) 



and is the medium 4- velocity. 

The result ([M]) describes the macroscopic evolution of 
a perturbative QGP in the presence of J'^ , which acts 
as a source of energy and momentum. In other words, 
J" gives the distribution of energy and momentum de- 
posited into a perturbative QGP due to the passage of 
a fast parton. This is precisely the quantity I set out to 
evaluate. In the remainder of this work, I will evaluate 
the source term both with and without dielectric 
screening and examine the resulting distribution. 
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III. EXPLICIT EVALUATION OF THE SOURCE 
TERM 



I begin by simplifying (|23p with an integration by parts 

(26) 

Po 



(2^) 



where 5qi, is the Kronecker-Delta symbol and, as men- 
tioned previously, j range from 1 to 3. I specify the 
medium to be a locally thermal plasma of masslcss gluons 
with the distribution 



/o 



2{Nl 



1 



S/3p" - 1 



(27) 



where Nc is the number of colors and = T is the local 
temperature of the medium. With this specification the 
distribution, /q, now contains the only dependence upon 
the magnitude of the momentum, p'^ ~ p, in (|26p . It 
follows that 



(27r 



Po 



47r 



((^oi/Vj + (5i,.)-D' 



(28) 



where v = p/p" is the directional velocity of medium 
particles. 

Inserting the explicit form of Dij gives 



r{x) 



^ — ik-x—ik' 



dt' I d^kd^k'dw^x 
An 



(29) 



where ttid = \/ Nc/3gT is equal to the Debye mass for 
gluons in the HTL approximation, and I have introduced 
the Fourier representation of the color Lorentz force 



d'^ke-'^-''F^{k). 



(30) 



The generalization to include Nf flavors of 
quarks/antiquarks in the medium can be done in a 
straightforward way by adding the appropriate thermal 
distribution to (P7)) 



2(iV2 - 1) ^ m^Nf 



(31) 



and repeating the steps which lead to ([29|) . The 
only modification is that the expression for the Debye 
mass appearing in (|29p now includes flavor: mo = 
\J {2Nc + Nf)/6 gT. In the numerical results that follow 
I have set Nf = 0. 

I next make the assumption that the typical momen- 
tum transfer between the source particle and the medium 
is small compared to the typical momentum of a medium 
particle. This assumption allows a straight line approxi- 
mation relating the position of a medium particle at time 
t' to its position at time t: 



(32) 



x' =x + v(t'-t). 



The t' integration can be evaluated after plugging p2| 
into (|29p and introducing an exponential damping factor 
ie, yielding 



(27r)8 



V ■ E°(fcO ((5o.v ■ E°(fc) + S,, {Efjk) + (v X B)f (fc))) 
47r(w' - k' • V -f k) 



(33) 



Next consider the integral over v. It can be evaluated 
by choosing a frame in which k' — k'z, performing the 
integral, and then rotating back into a frame in which 
k' is arbitrary. The explicit details of the calculation are 
given in Appendix B. The result is most easily expressed 
by dividing the source term into two pieces, and J2, 
where J'^ + J2 = J'^ ■ After integrating out the delta 
functions in ([l7]-[Tll) I find 

Jr(x) =^ J d3fcd3fc'eHk+k') (-uO^^(k')(^o.E'^(k') 
•E''(k)+5,,(E'^(k') X B''(k)),) 

(34) 



and 



(27r)» J 

X (E''(k') • k') (^(r2i(k')u • k' - (k' 
(5o.(E'^(k).k')+(5..(k' xB'^(k)), 



where 



k' + u-k' 



u k' 



2k''- 



(35) 



iTTU ■ k' 1 

—J^ - T7 (36) 
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1 



In 



k' + u- k' 



{k" - (u • k') 



l\2\ 



k' - U - k' 

- (u • k')^) + 2fc'u • k') 



(37) 



In the above expressions k' is the unit vector in the k' 
direction, k' is the magnitude of k', u is the source par- 
ton's velocity. The exponential damping factor in (|33p . 
le, sets the sign of the imaginary terms in (|36p . (j37p . 

At this point it is convenient to specify a direction 
for the source particle's velocity which I do by choosing 
u = M z. I also choose to work in plane polar coordinates, 
kr and 0, such that kx = kx cos (j) and ky = fc-r sin (/). 
With these choices the only dependence upon the mo- 
mentum space variable (j) is in the exponential and terms 
proportional to cos (j) or sin (j). The exponential depends 
upon (j) through the term ikx (x cos (f> + y sin 0) which can 
be rewritten as ikTpcos[(j)— a], where x = pcosa and 
y = psina. With these simplifications the entire (j) de- 
pendence of and can be integrated out by using 
the relations 



2tt 



2tt 



exp [ikxpicos [(f> — a])] 



1 

cos ( 
sin f 



Joipkr) ' 
fJiipkr) 



(38) 



where J,i(x) is a Besscl function of order n. The cal- 
culation is further simplified due to the fact that the k 
and k' dependence can be completely factorized, reducing 
the remaining four-dimensional integration to a product 
of two-dimensional integrations. Using the explicit field 
dependence listed above, I find there arc 12 unique two- 
dimensional integrations which must be performed. They 
are given by 



.^12. 



dkzdkr 



— I dkzdkx 



kDL{k) 



JoipkT)kTk^' 



dkzdkr- 



— / dkzdkr 



k^Drik) 

^i{z—ut)kz 

DL{k) 

^i{z—ut)k:, 

k^Drik) 



jQ{pkT)kTkl 

Jo{pkT)k^ 
Ji{pkT)kzkr 

■h{pkT)kzk^ 

■T)kl 

,2 

Jo{pkT)kzkT 

'Ji{pkT)k^k'^ 
Ji{pkT)k^ 



Ji{pkrr)k^k^ 
Ji{pkT)krp 



Ji{pkT)k'^kl 
Jo{pkT)kzk^, 



(39) 



where ^^(fc) ee k^eL{k), Drik) = k^ - u'^klerik). 
In terms of the above definitions the source term is 



g^T^iQi? x{y) 
167r4 p 

X (w^Cseg + (^1^7 



(40) 



'aai)) 



and 



167r4 



167r4 

u«6+e4)(e7 + w2ai)) 



(41) 



(42) 



The problem is now reduced to the evaluation of the 12 
a terms. Before considering the full results of such an 
evaluation it is useful to calculate the above expressions 
without including the medium's response to screen the 
fields. I perform this calculation in the next section and 
obtain an analytical result. 



IV. THE UNSCREENED SOURCE TERM 

The unscreened source term is found by setting the 
dielectric function given in (PP]) to unity, or, equivalently, 
setting DL{k) = fc^ and Drik) = k^ - u^kl in It is 
easy to verify that with this simplification and 
can be rewritten as 



J 



^^^^-^^^i|^^6(«^a + (a-«^6)), (43) 



-«2-^^'e5)(e8-^^'a2)) 



A quick perusal of the source term given by (|41l43l44p 
reveals that the only combinations of •■•j ^2 that need 
to be evaluated are (^ — ^^^12) and <^g. These relevant 
combinations are given by 



dkzdkx 



^i{z — ut)k~ 

kl + 72^2 



■y^Ji{pkT)k^ 




^9 


Jo{pkT)kzkT 




a - u^^i2 



(45) 

where 7^ = (1 — u^) ^. 

The above expressions can be evaluated in a straight- 
forward manner by first performing a contour integration 
in the k^ variable. The general form needed is 



dk. 



^±ik~ (z—ut) 



{kl + /fc2 



'I'kl 
kz 



ire 



^{z—ut)kT'y 



1 

±i 



(46) 



where ± refers to the sign of (z — ut). Similarly, one can 
use the relation 



dkrkre-^-^^^'-'''^ 
1 



Jiipkx) 
Joipkr) 



(/?2+72(z-ut)2)3/2 

to obtain the result 

7r7 



P 

7I2 — ut\ 



(47) 



a 

$8 - u'a 



(p2+^2(^_yi)2)3/2 



P 

i(z — ut) 



(48) 
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Before attempting an explicit evaluation of , . . . , it 
is worthwhile to consider the general form of the source 
term given by (|41|43l44p . It is interesting to note that 
(Plj) and ([33]) have a similar form. If one were to re- 
place ^1 — > (u^^2 — w^^s) in ((4T|) then one would find the 
relation 



J, -> uJq 



(49) 



I am motivated to find a relationship between , ^2 , and 
^5 that would allow me to write the source term in a form 
similar to that given by I can (trivially) write 

a = (^"6 - u^a) + (a - ^'6) + (50) 

which allows 

J = 4 [u a <6 + ii) - [[u a - u a) 

(51) 



and 



riv) ^ 2x(j/)^^^7 ("'a + (6 - u%)) (52) 



where the function d is given by 



d{p,z,t) 



87r(p2 4-^2(^_yi)2)3/2- 



(53) 



Notice that J^^J*) now has the same form as — uJ'^ 
apart from the factors of x and {z — ut). The source 
term can be written in the following compact form: 



where 



r{x) = (J"(x),uj°(x) - Jv) 



27 



(54) 



Jv = (x - ut) ^d{p, z, t) {u^^5 + (a - "'6)) • (55) 

I now continue by examining the terms ^7 ■■■iS.6 which 
are considerably more difhcult to evaluate than •■•j C12 
because of their dependence on the two Qi terms. The 
result and method of their evaluation is discussed in de- 
tail in Appendix C. Inserting the ^ into ([331) and ([55)1 I 
eventually find for the source term 



J°(p, z, t) =d{p, z, t)'-fu^ 

X I 1 - ^^7- ^ I 2_ + 



"fup 



{zl 4- p2) 



p2 + 2;2 ^2 



3^{p,z,t) =(x- ut) 



Sir {zl + p2) 

( U^p^ + (z2 ^2 ^ I 2z2 + 



,2^ (9.2 ^ ("'+2)P' 
^2 



V 



(^272 + ^2)^ 



2uz- 



(57) 

where z_ = (z — ut). 

Equations ([55)) and ([57)) . combined with ([M)) . give the 
distribution of energy and momentum deposited into a 
pcrturbativc QGP per unit time by a parton with con- 
stant velocity u = w z at position x — utz, in the absence 
of screening. That this result can be expressed in rela- 
tively simple closed form expressions is remarkable. The 
result for J°{p,z,t) is plotted in Fig. [Tjfor two different 
values of 7. The total energy deposited per unit time 
(or length) can be obtained by integrating ([56)) over all 
space: 



dE 
dx 



dxj"{p,z,t) 



In 



72^ 



(58) 



where y = cosh" "'^(7) is the rapidity and pmax and pmin 
are infrared and ultraviolet cut-offs for the p-integration. 
This result matches the standard leading-logarithmic 
result [13 for collisional energy loss in the relativis- 
tic (u —f 1) limit if one chooses pmax = I/'tt-d and 
Pmin = l/(2^£'pT), where Ep is the energy of the fast 
parton. It is interesting to plot the u dependence of ([55)) 
which is done in Fig. [2] The total momentum deposited 
per unit time can likewise be obtained by integrating 
3{p, z, t) over all space (recall that J = uJ° — Jv)- One 
can verify explicitly that the source term satisfies the 
energy-momentum relation 



J dxu ■ 3{p, z,t) ~ J dxj'^{p, z,t). 



(59) 



A careful observer may notice that the denominators of 
([56)1 and ([57)) contain terms involving both [z^'f^ + p^), 
which encodes the Lorentz contraction of the fast par- 
ton's field configuration, and i^z"^ + p^) , which is sensi- 
tive to the rest frame of the medium. One may ask why 
the result, in the absence of dielectric screening, should 
depend on the rest frame of the medium. The answer 
can be traced back to the t' integration in ([^^ which 
records the history of the fast parton's interactions with 
the medium. In that expression the fast parton inter- 
acts with a medium particle at t' , which is expressed by 
Fj{k'), and again at t, which is expressed by Fi{k). In 
between t' and t the medium particle propagates with a 
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for 7 = 20 




-2 




FIG. 1: (color online) Three-dimensional plots of the scaled energy density deposited per unit time by a fast parton moving 
in the positive z direction for as = I/tt and Nc = 3 [see 156). The plots are scaled by rn^{Qp)^ « 20.4 GeV/fm^ for a gluon 
moving in a QGP at 200 MeV. As one might expect, the distribution becomes Lorentz contracted for large values of 7. (In the 
above plots, I have set y — 0.4 GeV^^ to avoid plotting the origin where the source is divergent.) 



1- 



cosh ^(7) 



Atu^) with p = 1 [GeV"^] 



u y 




FIG. 2; Plot of the velocity dependence of —dE/dx as given 
by (|58[) where it is in units of c = 1. 



momentum that depends on the properties of the medium 
[see ([32]) ]. Thus the second interaction, at i, depends 
on both the field configuration of the fast parton and 
the properties of the medium, even in the absence of 
screening. One can observe this effect in Fig. |3l where 
J°(p, z, t)/{'-fu^) is plotted as a function of 72 for fixed p. 
If J°{p,z,t) depended only on the field configuration of 
the fast parton then the curves shown in Fig. [3] should 
be identical for any velocity. However, in the case of a 
highly relativistic source parton one sees a slight shift 
in the curve reflecting the dynamics of the medium in 
between the two interactions at t' and t. 



For a highly relativistic source parton Eqs. ([55)1 and 
57)) can be written in a simplified form by considering 



that terms with (2^7^ + p^) in the denominator will be 
severely damped except for z w ut. In this case, I keep 
(z — ut) to first order (and drop terms proportional to 




Y (z-ut) [GeV ] 



FIG. 3: (color online) A plot of [p, z,t) / {-yv?) as a function 
of 7(2 — ut) for the parameters u = 0.99, ttid = 0.4 GeV, 
Ua = I/tt and p = 1 GeV~^. If J°{p,z,t) depended only on 
the field configuration of the fast parton the curves should be 
identical. As one can see, there is a slight shift in the curve 
when going from 7 = 1.25 to 7 = 20. This shift reflects the 
dynamics of the medium in between the two interactions at 
t' and t (see discussion in the text). The fact that the shift 
in the curve is small even though 7 has changed by a factor 
of 16 suggests the effect is small. 



1/7, etc.) and the result becomes 



J°{p,z,t) = d{p,z,t)ju^ 1 



3v{p,z,t) = {x. — ut) d{p, z,t) 



7u z_ 



z^7^ + /O^ 



^17^ + p^ 



(60) 
(61) 



In the next section, I consider the numerical evaluation 
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FIG. 4: (color online) Plots of the scaled energy density deposited per unit time by a parton moving in the positive z direction 
with velocity u = 0.99 and Qs = I/tt, both with and without medium screening (compare with Fig. [l]). There is a very similar 
structure in both cases; however, as one would expect, the screened result shows some damping in the transverse direction. 




FIG. 5: (color online) Plots of the scaled transverse momentum density deposited per unit time both with and without medium 
screening included for the same parameters as in Fig. 2] The structure shows noticeable differences in the two cases. In 
particular, the screened result has a positive bump just behind the source parton which is absent in the unscreened result. 



p = 1.0 [Ge\rM 

Screened 

Unscreened — 



[GeV^] 







— — 

-3 -2 -1^ 

-o.ooa 


■ 


2 


3 


-0.004| 








- 0.006 









(z - u t) 



FIG. 6: (Color online) Plot of transverse momentum density 
deposited per unit time by a gluon both with and without 
medium screening included (compare with Fig. [S}. The plot 
is for u = 0.99, ttid = 0.4 GeV, and Qs = I/tt. 



of the source term with medium screening included and 
compare to the analytical result without screening. 

V. RESULT WITH SCREENING 

The evaluation of the source term (|40l41l42p with 
medium screening included must be done numerically. 
As mentioned before this amounts to evaluating the 12 

terms listed in . A discussion of the specific tech- 
niques used to perform the numerical integrations is given 
in Appendix D. 

In this section all evaluations are for a parton mov- 
ing along the positive z axis with velocity u = 0.99 and 
I have set mD = 0.4 GeV, A'c = 3 and = I/tt (the 
choice of w = 0.99 corresponds to 7 « 7, or « 3 GeV 
for the source gluon, which is likely too small a choice 
for the fast partons produced at RHIC but lends itself 
nicely to plotting the result). The result for the energy 
density deposited per unit time both with and without 
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medium screening is given in Fig. [4] One sees a very 
similar structure in both cases except that the screened 
result is moderately damped in the transverse direction. 
A similar plot is shown in Fig. Ofor the transverse mo- 
mentum density deposited per unit time. Here one finds 
more noticeable differences in the two cases. For instance, 
the screened result has a positive bump just behind the 
source parton which is absent in the unscreened result. 
This would correspond to matter being pushed outward, 
away from the source parton, in the vicinity of the bump. 
The differences are further highlighted in Fig. [HI Plots 
of the parallel momentum density distribution, J^, are 
not included because they are very similar to J° for a 
relativistic source. 



VI. DISCUSSION 

In this paper, I have derived the space-time distribu- 
tion of energy and momentum deposited by a fast parton 
traversing a perturbative quark-gluon plasma by includ- 
ing physics at the Debye screening scale. The approach 
has been to treat the fast parton as the source of an exter- 
nal color field perturbing the medium. The final result 
depends on two parameters: the source strength asQp 

and the Debye mass tod = y/Nc/3 gT. 

Several assumptions were made in deriving the result. 
The first, and primary, of these is that the medium is 
perturbative in the coupling as. This assumption is re- 
flected in the derivation of p5|) , the expansion of the path 
ordered gauge connection, Uabix, x'), and the choice of 
fields pn - [T9|) . Another assumption is that the medium 
is in local thermal equilibrium. This is the basic assump- 
tion of hydrodynamics and is necessary to obtain (|24|1 . I 
also assume that the typical momentum transfer between 
the source particle and the medium is small compared to 
the typical momentum of a medium particle. This as- 
sumption allows the straight line approximation given 
by ((32)) . To get an idea of where this last assumption 
breaks down one can use the equipartition theorem for 
masslcss particles to find a relation between the typical 
momentum of a medium particle and the temperature: 
p = ST. For T = 350 MeV the typical momentum of a 
medium particle is about 1 GeV which places an upper 
limit on typical momentum transfers. It follows that the 
straight line approximation breaks down somewhere on 
the order of fc > 1 GeV~^ (or distances less than 0.2 fm). 

In the absence of the medium's response to screen the 
fields, it was possible to obtain an analytical result which 
is given by (j56p and (|57p . This result shows many similar- 
ities to the result obtained when the medium's response 
to screen the color field is included. The most noticeable 
differences appear in the transverse momentum density 
distribution. The total energy deposited per unit time is 
obtained by integrating over all space and is 

given by ([55)1 for the unscreened result. 

The reader interested in reproducing the numerical re- 
sults for the screened source term should consult Ap- 



pendix D, where the specific techniques used to perform 
the numerical integrations are discussed. 
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Appendix A: Derivation of (|17l- 119^ 

In the HTL approximation, the source fields obey 
Maxwell's equations, which are given by 



V-D(x,<) = p(x,t) 

o»B(x, t) 



V X E(x,t) = -- 
V-B(x,0 = 



dt 



(A-1) 

(A-2) 
(A-3) 



VxB(x,i) = j(x,t) + ^5^Jl^ (A-4) 



dt 



where the source density is 



p(x, i) = gQpS{x-ut) 
j(x,i) = u/9(x,i) 



(A-5) 
(A-6) 



(in the above equations u is the source particle's velocity 
I have suppressed the color index, which I will restore in 
the final result). It is convenient to solve for the fields 
by taking the Fourier transform of Maxwell's equations. 
Using the general rule 

= (2^)^ / / ^''"""^"'^e^' ^) 
it is found that 

ik-D(k,cj) = p(k,Lj) (A-8) 

ikxE(k,tj) = iujB(k,uj) (A-9) 

ik-B(k,w) = (A-IG) 
zkxB(k,w) = j(k,u;) - iu;D(k,tj) (A-11) 

where 

p(k,w) = 2TrgQpS{uj-u-k) (A-12) 

j(k,c^) = up{k,Lo). (A-13) 

The function D is related to E through the dielectric 
tensor e,, : 



D, = e,,E, 



(A-14) 



(the summation over j is implied and I now suppress the 
(k, a;) notation). In an isotropic medium the dielectric 
tensor can be decomposed into transverse and longitudi- 
nal parts such that 



fc2 



fc2 



(A-15) 
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and D is now written as 

D = elEl + etEt 



(A-16) 



where El = k(k-E) and Er = E — El. One should note 
that k • Et = and k x El = 0. Maxwell's equations 
now take the form 

zcLk-EL = p (A-17) 

ik X Et = iiuB (A-18) 

ik-B = (A-19) 

ik X B = up - iLu{eLEL + etEt). (A-20) 

An equation for B can be obtained by crossing k into 
(|A-20P which combines with (|A-18p and (|A-19p to yield 
(with color index restored) 



^ _ 27rzgQg(k x u) 



B° = 



6{lu — u • k). 



One can then obtain Et by crossing k into (jA-lSP 



(A-21) 



E^ = — B X k 



2TrigQ^uj (uP - k(k • u)) 



Combining (|A-21|) and (|A-22|) with (|A-20|) yields 



6{uj — u • k). 
(A-22) 



i ( k^ 
E2 = (-1) up + i — Et - iujeT'Er) 



2^^5Q^(k • u) 



(A-23) 



ELiok'^ 



-'k5{uj — u • k). 



These results can be back transformed to position space 
to yield ([17] -[HI). 



ensure that only survives in (jB-2p . Making the re- 
placement dv ^ J dOdcj) sin 9 gives 



sm f cos ( 



47r([j' — k' cos 9 + ie) 



Similarly, in (|B-3[) the integral over d4> ensures only I^x, 
lyy, and Izz contribute: 



Ixx — lyy — . 



1 f' d^i-e) 



and 



1 {lj' - k'£_ + ie] 



^2 



The ri, terms are evaluated with the result 



k — uj 



k 



- 2 



(B-5) 



(B-6) 



(B-7) 



k — UJ 

k + iu 



k-i 



Equations (|B-4I - IB-6P can be rewritten as 

Ij — Sjz^l 



(B-i 



-Qiik). 
(B-9) 

(B-10) 



Appendix B: Explicit evaluation of the v integration 

In going from Eq. (|33p to and (|55|) I performed an 
integration over v by choosing a frame in which k' = k'z, 
performing the integral, and then rotating back into a 
frame in which k' is arbitrary. In this appendix the ex- 
plicit details involved will be given. The relevant starting 
point is 



V ■ E'" {Sq^v ■ E° + 4. {Ef + (v X B)f)) 
47r(u;' - k' • V + ie) 



(B-1) 



where E'° is short for E''(fc'), etc. There are two basic 
integrals to consider, which have the form 



47r([j' - k' • V + ie) 



4:Tr{uj' - k' • V + ie) 



(B-2) 
(B-3) 



Now specify that k' = k'z and consider that since 
V = (sin 9 cos cj), sin 9 sin (f), cos 9) the integral over d(j) will 



Ijm = (SjxSmx + SjxSmx) ^2 + SjzS^z^S 
= Sjm^2 + SjzS,nz{-^^l — ^2)- 



(B-11) 



In the above expressions the subscript z denotes the 
orientation of the k' vector. Generalizing to an arbitrary 
k' is done by making the replacement Sjz k'^ in (|B-10[) 
and (|B-11[) . Going back to (|B-ip (and suppressing the 
color index) I find 



V • E'° (,5o.v • E-^ + 4. (Ef + (v X B)f )) 



47r(tj' - k' • V + ie) 

{EjEilj + eimkljmEjBk) — 



60. ((E' • E)f72 + (k' • E')(k' • E)(-rii - 

5,, ((k' • ^')E^VLi + 

f (E' X B),f72 -I- (k' • E')(k' X B),(^r2i - f^a) 
V k 



(B-12) 



which is the result given in ([M]) and (|35p . 
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Appendix C: Evaluation of ^i, ...,^6 



In Sec. IIVI I put off the explicit evaluation of the terms 
^1, ...,.^6 and directly went to the results and ([57)1 . 
In this appendix, I will give a detailed derivation of the 
expressions: 



= / dkzdkT 



dkzdkT 



e^i^-ut)k,ni{k) 
kDL{k) 

k^Drik) 



'MpkT)kTk'^' 
Jn{pkT)kTkl 

Jo{pkT)k^ 
.Ji{pkT)kzk^. 

Jo{pkT)kzk^ 

,2 1,2 



Ji{pkT)k^k'^ 



(C-1) 



Here, I am interested in the unscreened source term 
which is found by setting DL{k) = k^ and Drik) = 



i^fc^. I start by changing to polar coordinates (i.e., 
k and 9, where kz = kcosO and fcy = fcsin^) in which 
case the two terms can be written as 



Vli{k) 



2k 

G2{e) 

Ak 



(C-2) 
(C-3) 



where 

Gi{0) = 2MCOs6'tanh"^[ucos6'] - inu cos 9 - 2 (C-4) 
and 

G2{9) = 2tanh-^[ucos6l](l-u2cos2 6l) 



— Z7r(l — cos^ 9) + 2ucos( 



(C-5) 



Defining A(0) as some generic function, I now write the 
general form of the k integration as 



= / dkd9 sin 9e 



ik{[z — ut) cos Q) 



Jo{pksui9) 
Ji [pk sin 9) 



which can be evaluated using the relations [28| : 

1 



(C-6) 



dke^'^yjoikb) = 

y^b^ -{y + ie)2 



1 / 

dke'^yji{kh) = 1 



^\P- - (y + ief 




where 6 > and e is a positive infinitesimal quantity. 
Defining a = — ut)/ p and using the above integrals 
allows 



= / dfcd6'sin6ie*'=((^-"*)™''*) 
,1 



Jo{pk sin 9) 
Ji{pksin9) 



A,(0) 



d9- 
P 

d9- 

P 



■\J sin^ 9 — {a cos O^ia)'^ 



1 



p-\J sin^ 9-~{a cos 0-\-ie)'^ 



MO) 



U{9) 

m 



HO)- 



(C-9) 



where I have set t 
terms look like 



for conciseness. Explicitly, the 



d9-G,{9) 
2p 

^^^G^W 



U{9) 
U{9) cos^ 9 
L{9) cos 9 sin 9 

U{9) cos 9 sin^ t 
L{9) cos2 6'sin^ 



(C-10) 



where ^3 is omitted since it will not appear in the final 
expressions. 

It is imperative to consider symmetry before trying to 
evaluate the various terms. For instance, in ^1, ...,^4 the 
imaginary part of U [9) and the real part of L{9) only need 
to multiply against the term —iiru cos 9 due to symmetry 
considerations. Similarly, the real part of U{9) and the 
imaginary part of L(9) do not need to multiply against 
~iTTucos9. It turns out that due to symmetry the terms 
involving the imaginary part of U(9) and the real part of 
L{9) are easier to evaluate than those involving the real 
part of U{9) and the imaginary part of L{9). It makes 
sense to first get these simpler expressions out of the way 
and then focus on the more tedious ones. I introduce 
the notation ■Ci = ?ia + Cifej where now denotes terms 
involving the imaginary part of U{9) and the real part of 
L{9). Using Mathematica 6.0 these terms can be evalu- 
ated and I find 



^5 a 
^6 a 



7rau(3+a ) 

i-rvu / 1 _ z a(3+a 
' p X* 

_ z a(3+a^) 
P X'' 



(C-11) 



6p 



where x = Vl + a^. 

Now I turn to the more involved terms which in- 
volve the real part of U{9) and the imaginary part of 
L{9). It is convenient to note that 



Re[U{9)] = sin 9 Re 



1 



\/ sin^ 9 — cos^ 
1 

a/ sin^ 9 — a? cos^ 1 



(C-12) 
(C-13) 



For the sake of book keeping, I note there are 7 distinct 
integrations which need to be performed. These terms 
are given by the real parts of 



Vt= I d9 



smf 



\/ siv? 9 — cos^ 



(C-14) 



where 



























T7. 





1 

0082 6* 

COS 9 tanh^ ^ [u cos 9] 
cos'^ 9 tanh^ [u cos 9] 
cos 9 sin^ 9 tanh~ ^ [u cos 9] 
cos2 6'sin2 6'/ (1-^2 cos2 9) 
cos'' 9/{l-u^ cos^ 61) 



In terms of these the are 









1 








P 


-$66- 





Uip3 - (^1 
U(^4 - (^2 

((y55 + uipe)/2 



(C-15) 



(C-16) 



Using Mathematica 6.0 1 find for the real parts of the ip 
terms: 



(C-17) 



= — J ( X - + 1/7^ 



(X^ + 2^2) 7i_ 



7%4 \ 



TT / 1 1 
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(C-19) 

(C-20) 
(C-21) 

(C-22) 



(C-23) 



tP2 



2x3 



(C-18) Combining the above results I find 
I 



^5 



-TT ^-^/a^ + 1/7^ — au^ 



TT (2^a^ + ^^72 (a3«3 + x^) + (1 + a^^^) (2u2 + x2(l ~ Su^)) - Sx^) 



(2 + + x^ (x - 3a^^) + 2ya2 + ^(x^ - u')^ 



u6 - W ^5 



u^e + 2^4 



TTM ^^a (^-207^ + 2^72 y'a^ + + a) - 1 



2x%/a2 + ^72 



(^-a2 + (a^ + 3) ua^a"^ + + 2u'^ - 1^ 
4 /„2 I 1 



ITTU 



(C-24) 



which are combined with ([48|) to yield ([56|) and ([57]) . grations. The terms are given by 



Appendix D: Numerically Integrating Equation (|39p 

As mentioned in Sec. |Vl the evaluation of the 12 
terms listed in (|39p with medium screening included must 
be done numerically. In this appendix I will discuss the 
specific techniques used to perform these numerical inte- 



— \ dkzdkT 



kDL{k) 



Jo{pkT)kTkl 

JoipkT)k^ 
.Ji{pkT)kzkj 



(D-1) 
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dkzdkx 



JoipkT)kzk^ 
Ji{pkT)k?^kl 



— I dkzdkx- 



DL{k) 



dkzdkx 



^i{z—ut)k- 

k^Drik) 



Ji{pkT)kl ' 
Jo{pkT)kzkT 



Ji{pkT)k'^k'^' 

Ji{pkT)k^ 
Ji{pkT)k'^k'l 
.JoipkT)kzk^, 



(D-2) 
(D-3) 

(D-4) 



where the k^ integration runs from ±00 and the kx in- 
tegration ranges fr'om — > 00. Since each S,i contains a 
double integration which ranges over an infinite interval, 
care must be taken to cast each term in a form that will 
tend to zero as quickly as possible in order to make nu- 
merical evaluation feasible. In the case of the integrals 
over dkz , we can take advantage by the exponential term 
gi(z-iit)fej |-jy Jjgriding the contour in the imaginary plane. 
The strategy works as follows: start with a function of 
the form 



-^*^'''f{kz)dkz 



(D-5) 



where f{kz) is an analytic function. Next, bend the con- 
tour of the integral in the complex plane by making the 
variable change 



and 



dk, dre'"S'sn[r] 



(D-6) 



(D-7) 



where the limits of integration on r are from —00 to 00. 
Under this change, (|D-5P becomes 

cxp [i{z - wt)re''"S's"M]/(re'"S's"M)e*"'^'snW^^ 

_ J 

(D-8) 

and now the exponential term is damped in (|D-8p [pro- 
vided a has the same sign as (z^ut)] making the integral 
more numerically manageable. (|D-8[) and (|D-5p will be 
equivalent if the integrand vanishes as r — > ±00 and no 
poles from f{kz) are crossed when bending the contour 
(this places a restriction on the value of a). 

Before applying this technique it is necessary to con- 
sider where the have poles. Any poles located directly 
on the imaginary axis will not affect the result so we 
can ignore those. That leaves us to consider where the 
zeros of Dxiukz) and DL{ukz) are. An analytic solu- 
tion for the zeros of Dxiukz) and D^iukz) is impossible 
because of the logarithmic terms. However, the zeros 
can be searched for numerically. Using the parameters 
u — 0.99z and ttt-d = 0.4 GcV, it is found that all poles 
are sufficiently close to the imaginary axis so that the 
value of a used is not restricted for any reasonable value 



for terms with Driukz) in the denominator. I choose to 
use |a| = 1. However, for terms with DL{ukz) in the de- 
nominator, it is found that for values of fc^ ^ 1/2 poles 
become increasingly close to the real axis. For this rea- 
son, I choose to break up the integration in the terms 
with DL{ukz) in the denominator into a term 



dkx 



dk. 



(D-9) 



which I bend in the complex plane with a value of |a| = 1 
and a term 

dkr / dkz (D-10) 

J —oc 

which I do not bend in the complex plane. For the terms 
I do not bend in the complex plane it is usually necessary 
to subtract off the asymptotic form of the kz integration 
in order to make the numerical integration efficient. The 
asymptotic form must then be added back. Similarly, for 
the terms I do bend in the complex plane, I usually must 
subtract off the asymptotic form of the fc^ integration to 
provide efficient convergence. 

In the case of z = ut, it does not make sense to bend 
the kz integral in the complex plane since there will be 
no damping term from the exponential. In this case, I 
instead perform the numerical integration by subtracting 
off the unscreened form of in the integrand (which has 
the same large k form) and then re-adding after evalua- 
tion. 

To make these ideas more concrete, consider ^7 for the 
case oi z ^ ut. I write 



dk. 



dkre'^''-'"'^^' Jiipkr) 



1 



DL{ukz 

(D-11) 

where I have subtracted 1 to make the kx integration 
more efficient. I bend the above contour in the complex 
plane and perform the integration and then re-add 

dkz I dfcTe'(^-"*)'=-Ji(pfcT)-(5(z-Mi) (D-12) 
-00 Ji 

which can be ignored, since it only contributes when z = 
ut. Next, I integrate 

dkz C dkre'^'-^'^'^' Ji{pkT) ^ f, , 
-1 Jo DL[ukz) 

1 / 7,2 



(D-13) 

where I have subtracted k^ / k1 to make the kz integration 
more efficient. I then must re-add 
»i 



2Re / dkz / dkre'^^-""'^^^ .h{pkT)-^ 

Jl Jo . 

Up) 



(D-14) 



(2(cos(z)-t-zSi(z))-7r|z|). 
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Considering ^7 for the case oi z = ut, I write 



and evaluate at z == ut. 



/oc poo 
-00 Jo 



(D-15) 

and then analytically re-add [see (|48)) ] 



'^T + «3 (D-16) 







(p2 + (z _ Mt)2)3/2 
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